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Abstract 

We provide a self-contained theoretical analysis of the dynamical response 
of a one dimensional electron system, as confined in a semiconductor quan- 
tum wire, within the random phase approximation. We carry out a detailed 
comparison with the corresponding two and three dimensional situations, and 
discuss the peculiarities arising in the one dimensional linear response from 
the non-existence of low energy single-particle excitations and from the linear 
nature of the long wavelength plasmon mode. We provide a critical discussion 
of the analytic properties of the complex dielectric function in the complex 
frequency plane. We investigate the zeros of the complex dielectric function, 
and calculate the plasmon dispersion, damping, and plasmon spectral weight 
in one dimension. We consider finite temperature and impurity scattering 

effects on one dimensional plasmon dispersion and damping. 
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I. INTRODUCTION 



In an important recent experiment, Goiii et al. [l[] measured the wave vector disper- 
sion of the collective charge density excitation ("plasmons") in a one dimensional electron 
gas confined in GaAs quantum wire structures by using the inelastic resonant light scat- 
tering spectroscopy. The measured one dimensional plasmon dispersion agrees remarkably 
well with the earlier calculated random-phase-approximation(RPA) predictions 0, which 
initially was considered surprising in view of the widespread belief that one dimensional 
electron systems are singular Tomonaga-Luttinger liquids where standard Fermi liquid type 
electron gas theories (of which RPA is a prototypical example) are inapplicable. The quan- 
titative agreement between RPA predictions and experimental results was later explained 
by the fact that the calculated RPA plasmon dispersion and the Tomonaga-Lutinger the- 
ory for the collective charge density excitations (Tomonaga-Lutinger bosons |4[]) of the one 
dimensional electron system are equivalent by virtue of the vanishing of vertex corrections 
H in the irreducible polarizability of the one dimensional electron gas. Subsequent to the 
original experimental discovery of the one dimensional plasmon by Goni et al. |l| several 
theoretical papers have appeared in the literature explaining various special features of the 
experimental data 0-0. 

The purpose of this paper is to provide a complete, self-contained theoretical analysis of 
the dynamical response of a one dimensional (ID) electron system based mostly on RPA. 
It is surprising that in spite of considerable interest in the subject matter such a compre- 
hensive analysis does not exist in the literature. By providing rather complete results for a 
well-defined simple model and by extensively comparing the one dimensional results to the 
well-studied two and three dimensional RPA results, we hope to dispel some misconceptions 
and erroneous claims made in the recent theoretical literature and to accentuate the pecu- 
liar aspects of one dimensional plasmons (as against the corresponding higher dimensional 
situations) arising from the severe phase-space restrictions in one dimension. We emphasize 
that RPA takes on particular significance in one dimension because its domain of validity 
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for plasmon dispersion is large (up to second-order in the wave vector in one dimension as 
against only first order, as implied trivially by the f-sum rule, in two and three dimensions). 

The rest of the paper is organized as follows: In Sec. II we provide the basic RPA 
formalism and applicable formulae; in Sec. Ill we give the numerical results for plasmon 
dispersion and damping, and in Sec. IV we provide the dynamical structure factor (which 
is a direct measure of the light scattering spectral intensity); in Sec. V we go beyond RPA 
by including some local field corrections in the Hubbard approximation, and we conclude in 
Sec. VI providing a summary of our results. 



II. THEORY 

The elementary excitation spectrum of an electron gas is given by the dynamical structure 
factor S(q, uj) which is proportional to Im[e(g, uj)^ 1 ] where e(q, uj) is the dynamical dielectric 
function of the electron system. Within RPA, the dielectric function e(q,oj) is given by 

e(q,u) = l-v c (q)Ho(q,u), (1) 

where v c (q) is the Coulomb interaction (in wave- vector space), and U (q, uj) is the leading- 
order irreducible electron gas polarizability function {i.e. the so-called bare bubble or the 
Lindhard function in the relevant dimension, Fig. 1(a)). In an ideal ID electron system the 
Coulomb interaction v c (q) in the wave vector space is logarithmically divergent. By using a 
more realistic finite-width ID model flJlOll, the Fourier transformed ID Coulomb interaction 



is given by 

v(q,y-y') = —K (q\y-y'\) (2) 

where e& is the background lattice dielectric constant, and K is the zeroth-order modified 
Bessel function of the second kind. In the finite-width ID model we use a square well 
potential with infinite barriers at y — —a/2 to y = a/2. The ID electron gas is along the 
x-axis and the wire width in the y (z) direction is taken to be a (zero). In typical quantum 
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wires, the wire thickness along the growth direction is much smaller than that in the lateral 
direction, justifying our ID finite width model. The Coulomb matrix element in this square 
well potential becomes 



v c (q) = —[ a dy f dy'K (q\y-y'\)cos 2 (^-)cos 2 (^-). (3) 
Its asymptotic forms [11 
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for \qa\ — > oo, 



v c {q) = { e T ( 4 ) 
^[K (\qa\) + 1.9726917--], for \qa\ -»• 0, 

where a is the cross-sectional confinement size of the quantum wire. The noninteracting 
polarizability can be exactly calculated in any dimension, and we show here the results 
for the one dimensional case only. (The corresponding two [0 and three |y| dimensional 
results can be found in the literature.) The ID polarizability function n (g, lj) is 

n o(g , u) = (2s + 1) J c ^ E{k + q) _ m _ u , (5) 

where s = 1/2 is the spin of the carrier, E(k) = k 2 /2m — // with chemical potential //, and 
f e q(k) is the Fermi distribution function with k as the ID wave vector. (Throughout this 
paper h = k B = 1.) 

Plasmon dispersion and its lifetime or damping is determined by the zeros of the complex 
dielectric function such that 

e(q,z p )=0, (6) 

where z p is the complex frequency (z p = u p — ia), where u p and a are taken to be positive 
functions of the ID mode wave vector q. (This ensures stability of the system.) The condition 
Re[,2p] = uj p (q) gives the plasmon dispersion and Im[z p ] = a(q) gives the damping rate of the 
plasmon mode. Since a is positive, the analytic continuation of the dielectric function into 
the lower frequency half plane is needed. Then the polarizability for the complex frequency 
[z = us — ia) becomes 
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Without any damping (i.e. when a = 0) we find the dispersion relation for an undamped 
quantum wire plasmon at T = to be 

- A(q)ul(q) ~ u 2 _(qV 1/2 



u p (q) 



A(q) - 1 

where A(q) = exp[q-K /mv c (q)]. In the long wave length limit (q — > 0) this mode becomes 
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uj p (q) = q 



v 2 F + -v F v c (q) 



71 



1/2 



0(q 3 ). 



(11) 



This RPA plasmon dispersion up to second order in q is the same as the dispersion of the 
collective charge density excitation in the exactly solvable Tomonaga-Luttinger ID model 



Since u) p (q) (Eq. (10)) is greater than tf + (g), which is the upper boundary of the electron- 
hole pair continuum for all q, we expect this mode not to decay to electron-hole pairs within 
RPA. In three dimensions, however, u; + (q) becomes larger than the plasmon energy at a crit- 
ical wave vector, q c m u p /vp. If q > q c , e(q, u) has a finite imaginary part and the plasmon 
mode enters the pair continuum regime, becoming damped consequently (Landau damp- 
ing). The plasmon mode inside the Landau damping region decays by emitting single pair 
excitations which is now allowed by energy-momentum conservation. For small Im[e(g, u)], 
this damping rate can be calculated by taking z = uo(q) — ia in Eq. (1) (with u as a real 
frequency, and a/u <^ 1) |0|, to obtain 



Re{e(q,u p (q))] = 1 - v c (q)Re[U (q, LO p {q))\ = 0, 



(12) 
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Eq. (12) determines the dispersion relation u) p (q) of the plasmon, and Eq. (13) gives an 
explicit formula for the plasmon damping rate. In 2D, the plasmon mode approaches the 
pair continuum and meets the damping region at some critical wave vector, q c . Its difference 
from 3D is that the plasma mode simply ceases to exist at the upper boundary of the pair 
continuum at q c . Since the plasmon does not exist inside the pair continuum, the analytic 
continuation of the dispersion relation, which is useful to explain the damping of the 3D 
mode, is invalid in 2D. Within RPA the 2D plasmon mode does not decay by single pair 
emission as long as it exists. It simply disappears at q c . A similar phenomenon is even more 
striking in ID. Since the plasmon mode always lies above the upper bound of the ID pair 
continuum, the ID plasmon is undamped within RPA in a pure system. Because of this 
fundamentally undamped character of ID RPA plasmons, there are no solutions of Eq. (6) 
with a finite a in ID. In some sense one can define a "true' collective mode by the zeros 
of the real part of the dielectric function, Re[e(g, u)} = 0. Within RPA there are always 
two zeros of Re[e(g, u)\ for all (ID, 2D, 3D) dimensions - one inside the pair continuum 
and the other outside the pair continuum. Thus, one may be misled into concluding that 
there are two collective modes corresponding to the two solutions of Re[e(g, to)] = 0. The 
damping of the mode existing inside the pair continuum is, however, very large, making it 
an experimentally unobservable and theoretically uninteresting entity. (Within RPA, the 
strength of the mode outside the pair continuum is typically a thousand times greater than 
that of the inside mode.) There has been some recent confusion in the literature || regarding 
this point, and incorrect claims of there being two collective modes (corresponding to the 
two solutions of Re[e(q, u>)] = as discussed above) in ID systems have been made. We 
emphasize that the two zeros of Re[e(g, u)}, one inside and one outside the electron-hole pair 
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continuum, exist in all dimensions, and only the zero lying outside the continuum is a true 
collective mode carrying any significant spectral weight. 

The expression for the ID polarizability Ho(q, oS) given in Eq. (7) assumes that the system 
is pure and free from any defects or impurities. In general, impurity scattering causes the 
electrons to diffuse instead of moving ballistically. Impurity scattering effects are usually 
introduced by including impurity ladder diagrams in the polarizability consistently with self- 
energy corrections in the electron Green's function (Fig 1(b) and (c)). With the impurity self- 
energy correction the dressed Green's function becomes (within the Born approximation), 

G(q, w) = ^7^1 7T, (14) 

UJ - E(q) + ^sgn(^) 

where r is the impurity induced elastic scattering time in the Boltzmann approximation. 
Then the expression for polarizability II 7 (g, to) including impurity scattering is given by 

n 7 (^)=/^ U [ E ^ U) , (is) 

7W ' ; J 2vril -v 2 II(£,g,u;)' V ; 

where v is the impurity disorder potential which is assumed to be short-ranged (isotropic) 
with its strength given by 
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(16) 
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where Np is the density of states per spin at the Fermi level, which is given by 2m/ (nkp) in 
ID. In this calculation we use a simple random Gaussian white-noise model for the impurity 
disorder potential 



v(x)v(x r ) = v 8(x — x). (17) 

In Eq. (15), the leading order polarizability kernel in the presence of impurity scattering, 
H(E, q, lj), is given by 

Tl(E,q,uj) = J2G(k 1 E)G(k + q,E + uo). (18) 



In doing the k integral we approximate JL4 



£[ ] — Np j dek [ ]. (19) 

k 

By using this approximation Eq. (18) becomes for q <C kp, < u <C £j? 5 and in the 
weak-scattering limit (Epr ^> 1) 



n (£ , 5 ,.) = 2.V ( - B(B + w)) 



1 + r^ + °(« 4 ' 



(20) 



u + i/r 

where D = vpr is the one-dimensional diffusion constant, and the Heaviside 9 function 
ensures that the two poles are in different half-planes. Finally, the polarizability ILy(q,u) is 

Da 2 It 

n 7 fe w) = N F - ,_ r 9 / ; T . (21) 



t D<7 2 /t 
t (o;+i/r) 2 



(w + i/r) 

This expression approximately describes the long wavelength and low frequency electronic 
density response for diffusive (Epr ^> 1) electrons. For very large values of r, Eq. (21) 
reduces to the correct expression for the polarizability, H (q,u), of a noninteracting ID 
electron system, 

2 2 

U (q,u) = N F ^, (for q « cu/v F ). (22) 

For finite r, in particular, in the limit q, to — > 0, IT 7 (g, cj) shows the pole structure charac- 
teristic of the diffusive regime, 

Note that Eq. (21) is only a long wavelength and low frequency approximation to the full 
polarizability, Eq. (15), which is, in general, unknown within this diagrammatic approach, 
and can only be determined by solving the extremely difficult integral equation defined by 
Eq. (15). Since the exact expression for Hj(q, u) within this diagrammatic approach is 
complicated, we use in the numerical calculation a particle-conserving approximation to Eq. 
(15) for arbitrary values of q and u, given by Mermin |DJ. In this relaxation time conserving 
approximation |TJ|], with 7 = 1/r, the polarizability is given by 

TT ( \- (^ + ^7)n (g,^ + i7) m x 

Aq,u}) ^ + * 7 [no(g,^ + * 7 )/n (g,o)r [M) 
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This approximate formula gives the same diffusive long wavelength and low frequency be- 
havior (i.e. Eq. (21)) as one gets from the diagrammatic approach. We have used this 
formula (i.e. Eq. (24)) for our full numerical calculation for all q and u, taking the impurity 
scattering induced broadening 7 as a constant phenomenological parameter. (Note that 
the plasmon broadening ot(q) is, in general, determined by, but different from the impurity 
broadening 7.) 

We next examine finite temperature effects on ID dynamical response. Unlike 3D metals, 
where the Fermi energy (Ep ~ 10 K) is much larger than experimental temperatures, mak- 
ing the zero temperature formalism very accurate even at the room temperature, artificial 
semiconductor structures of reduced dimensionality have Fermi energies which are typically 
comparable to experimental temperatures (e.g., Ep = 4AmeV ~ 50K), and therefore finite 
temperature effects should be included in the theory. This is because typically these sys- 
tems have much lower effective electron density and effective mass, and much larger lattice 
dielectric constant than metals. Finite temperature effects are included in our theory by 
calculating Tlo(q, 00) at T 7^ 0, using the trick introduced by Maldague JTB| for a 2D electron 
gas. Its comparative ID expression is given by JTT] 



U (q,uj;T,ii) = - — — ^n (g, u; T = 0, fi + Tln(x)) 

Je-v/ T (x + \y 

+ r— T ^n (g,^;T = 0,/i-Tln(a;)) (25) 

Jo (x + iy 

where fi = fi(n, T) is the finite temperature chemical potential to be obtained from the total 
number of particles n = J dED(E)f(E), using the ID electron density of states, D(E), and 
the Fermi distribution function, f(E). 



III. NUMERICAL RESULTS FOR PLASMON DISPERSION 

In this section we present our numerical results for ID plasmon dispersion and damping, 
comparing with the corresponding 2D and 3D results as appropriate. In Fig. 2 we show 
the plasmon dispersion and the Landau damping region as a function of the wave vector for 
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r, = 1.4. Here r s = (-r-) 1 ^ 3 — , (— W 2 — , and -§-— is the dimensionless electron gas density 
parameter |T3"1 , |T7|1 for 3D, 2D, and ID, respectively, with n and ao being respectively the 
electron density and the effective Bohr radius. For a GaAs ID system, r s = 1.4 corresponds 
to a density of n — 0.56 x 10 6 cm _1 , = 12.7, the effective mass of m = 0.067m e , and a 
Fermi energy of Ep = 4AmeV. We choose the quantum confinement size a = lOOA in all 
our results. While we have chosen typical values of n and a for GaAs quantum wires in 
the extreme ID limit, our qualitative results are completely independent of our choice of 
parameter values. In the long wavelength limit the plasmon dispersion has different behavior 
in different dimensions. In one and two dimensions the plasma frequency goes to zero as 



gy|ln(g)| (ID) and ^Jq (2D) as q — > 0, but in three dimensions the plasma frequency is 
finite in the long wavelength limit. The area bounded by the solid lines (the upper one, 
uj + = q 2 + 2q and the lower one, U- = q 2 — 2g, with kp as the wave vector unit) indicates the 
ID electron-hole pair continuum region. Unlike 2D and 3D, where the pair continua have 
only upper boundary, u + , the ID pair continuum has a low energy gap between and co>_. 
This gap arises from severe phase space restrictions imposed JTT| by ID energy-momentum 
conservation, and leads to the non-existence of low energy single particle excitations in ID 
Fermi systems. The dispersion of the ID plasmon mode lies above the upper boundary 
uj + of the pair continuum for all q. This means that the plasmon mode does not decay 
into electron-hole pair excitations and is undamped in a pure ID system within RPA. Fig. 
3(a) shows (in the absence of impurity scattering, 7 = 0) the zeros of the Re[e(q, u)] and 
Im[e(g, a;)] for q = 0.3k F in the (u — a) complex frequency [z = u — iot) plane using Eq. (7), 
which is an analytic continuation into the lower half energy plane. There are two branches 
(Fig. 3(a)) of zeros in Re[e(g, uj)\ for finite a up to a max = 0.38, whereas Im[e(g, uS)\ = 
along the a = line and for < uj < uj- and u+ < u < 00. The condition for the existence 
of a plasmon is Re[e] = Im[e] = 0, i.e. e(q,cj) = 0, and therefore only the point P (a = 
and uj = oj p (q)) in Fig. 3(a) corresponds to a true plasmon mode. This mode (the point P) 
is undamped with a = 0. For finite a, even if Im[e(g, a;)] changes sign (as it does along the 
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dotted lines) it does not pass through a zero because its change is discontinuous (that is, 
| line | > 0). Since the dynamical structure factor is proportional to Im[e] _1 when Re[e] = 0, 
we expect very small spectral weight in the region where Re[e] = and Im[e] 7^ (along the 
solid line in Fig. 3(a)). The structure factor is broad and incoherent, rather than a 5-function 
collective mode peak, along the solid line in Fig. 3(a) where Re[e] = 0, Im[e] = (See section 
IV) . By investigating the analytic continuation of the dielectric function into the lower half 
of the complex energy plane Bonitz et. al. |§ recently concluded that there is a collective 
mode in ID with a 7^ 0. We believe, for reasons explained above, that this conclusion is 
erroneous and ID dynamical response is qualitatively no different from 2D and 3D in this 
respect - this extra "mode" along the solid line in Fig. 3(a) with Re[e] = 0, Im[e] 7^ 
is overdamped and is not a collective mode in any sense. The collective plasmon mode 
corresponds to the zero of the total dielectric function (i.e. both Re[e] = and Im[e] = 0) 
and has zero damping (a = 0) within RPA in ID, and there are no additional plasma modes 
for finite a in a single component ID system. The experimentally observed additional mode 
of ref. [TJ is successfully explained quantitatively by introducing two-subband occupancy 
into the model as has recently been done in the literature ||. The additional mode arises 
from the interaction between two different ID subband plasma modes - one is the in-phase 
mode and the other the out-of-phase mode; a strictly one subband model including analytic 
continuation can not explain the existence of the additional mode. 

The inclusion of impurity scattering (7 7^ 0) causes collisional damping of the plasmon 
mode introducing level broadening in addition to any Landau damping arising from the 
decay of the plasmon to electron-hole pair excitations. This corresponds to complex zeros of 
the complex dielectric function with finite a. Fig. 3(b) shows Re[e(g, w)] = lm[e(q,u)] = 
lines in the complex frequency plane for q = 0.3kp and an impurity scattering induced 
broadening 7 = O.lEp. Compared with Fig. 3(a) the point P has shifted to finite a in 
Fig. 3(b), and consequently the ID plasmon has finite damping in the presence of impurity 
scattering. In Fig. 4 the plasmon dispersion (u>(q) = cu p (q)—ia(q)) calculated in the presence 
of finite impurity scattering (Eq. (24)) is shown for various impurity scattering rates, 7. In 
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Fig. 4 the curves with uj > give the plasma frequency or the real part (u p (q)), and those 
with uj < give the plasma damping or the imaginary part (a) of the complex zero solutions 
(a> (<?)), i.e. e(q,u(q)) = 0. From the figure it is clear that the plasmon is overdamped below 
a critical wave vector q < q c . The critical wave vector, q c , below which the plasmon does 
not exist due to impurity scattering effects, depends on the density and the scattering rate 

qc = \hk^Y\ (26) 

where K = ^Jmeo/8ne 2 . 

It is known [Tj| that the 2D plasmon is also overdamped in the long wavelength limit 



in the presence of impurity scattering. The plasma dispersion of the 2D system for small q 
and uj in the presence of a finite collisional broadening 7 is easily calculated on the basis of 



Eq. (24) (using the corresponding 2D noninteracting polarizability IIo ||12||): 

uj = + ^\/ —f 2 + SV2r s q. (27) 
The critical wave vector, q c , at which uj becomes completely imaginary in 2D is 

*=iA;- (28) 

We show the 2D plasmon dispersion for different value of collisional broadening in Fig. 5(a). 
The 2D results in Fig. 5(a) are qualitatively similar to the ID results in Fig. 4. As in ID, the 
plasmon damping a(q) decreases with increasing 7 in the overdamped region because of the 
suppression of the electron diffusion rate for larger 7. When the impurity scattering rate 7 
is small the 2D plasmon dispersion is hardly affected by impurity scattering. A quantitative 
comparison of Fig. 4 and 5(a) show that for a given range of collisional broadening 7, the 
plasmon dispersion is more strongly affected in ID than in 2D by impurity scattering effects. 
In 3D, the plasma frequency has a long wavelength gap, with the long wavelength plasma 
frequency being finite in contrast to ID and 2D systems. Impurity scattering, therefore, 
has no qualitative effect on the 3D plasma dispersion (Fig. 5(b)), only small quantitative 
corrections. For small 7, the 3D plasmon is not overdamped for small q, and, therefore, 
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there is no long wavelength critical wave vector, q c , in 3D below which the 3D plasmon 
vanishes. (On the other hand, as discussed earlier, there is a Landau damping induced large 
wave vector cut off, q c , for 3D plasmons, and the 3D plasmon ceases to be a well defined 
mode for q > q c , with q c being the wave vector at which the 3D plasmon dispersion enters 
the electron- hole single particle continuum - see Fig. 2, for example.) 

Fig. 6 shows the finite temperature effect on the plasmon dispersion. Unlike impu- 
rity scattering effects the plasmon dispersion is only slightly affected by finite temperature 
effects. Also, in contrast to impurity scattering effect where increasing 7 decreases the plas- 
mon frequency, increasing temperature actually enhances the ID plasmon dispersion. This 
upward shift in energy with increasing temperature is a well-known phenomenon in plasma 
physics |19| , and happens in all dimensions. Thus, effects of finite collisional broadening and 
temperature are qualitatively different on plasma dispersion. In Fig. 6(b) we show our cal- 
culated ID plasmon dispersion and damping in the presence of both non-zero temperature 
and impurity scattering. 



IV. DYNAMICAL STRUCTURE FACTOR 



The dielectric function e(q,u) is related |13| , |T7H to the dynamical structure factor S(q,u) 



by 



S(q,u) 



nv c (q) 



Im 



(29) 



where n is the ID density and v c is the Coulomb interaction, as in Eq. (1). S(q,u) is, 
therefore, the spectral weight of the dielectric function. The dynamical structure factor 
gives a direct measure of the density fluctuation spectrum of the electron gas. Thus, many 
experiments such as inelastic electron and Raman scattering spectroscopies directly measure 
S(q,u). Equivalently, S(q,u) is a measure of the spectral strength of various elementary 
excitations of the electron gas. An undamped plasmon shows up as a 5-function peak in 
S(q,u) indicating the existence of a simple zero of e(q,u). The damped plasmon (a ^ 0), 
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however, corresponds at best to a broadened peak in S(q, uj) - for larger broadening, as 
for q < q c in ID and 2D with impurity scattering, the plasmon is overdamped and there 
is no peak in S(q,ui). Note that S(q,u) oc Im[e(g, u;) -1 ], and for our purpose we mostly 
show results for Im[e(g, u;) -1 ] because v c (q) in ID depends on the details of the electron 
confinement model (cf. Eq. (3)). 

In Fig. 7 the calculated ID RPA Im[e(g, u;)" 1 ] is plotted for different wave vectors as a 
function of the frequency for a clean system at T = 0. When both Re[e] and Im[e] become 
zero (i.e. e(q, u) = 0, which defines the plasmon mode), the imaginary part of the inverse 
dielectric function Im[e(g, cj) -1 ] is a 5-function with the strength 



where w p (q) is the plasmon dispersion for a clean electron gas. The sharp 5-function like 
peak in Fig. 7 corresponds to the well-defined ID plasmon mode. This figure also shows 
the electron-hole single particle excitation (SPE) continua, which show up as weak broad 
(incoherent) structures in S(q,u). The continua in Fig. 7 have been enhanced by a factor 
of 10 to make them observable on the scale of the plasmons. Note that there is only one 
observed peak, the plasmon mode, in Fig. 7 with any appreciable spectral weight. In 
contrast to 2D and 3D systems, this remains true in ID systems up to rather high wave 
vectors (> 2kp) because of the severe suppression of SPE. This suppression of SPE continua 
(up to rather high wave vectors) is a specific characteristic of the ID electron system. (Note 
that there is only a single plasmon peak in Fig. 7 in contrast to the erroneous claims of 
ref. 0.) The low wave vector suppression of SPE is a generic feature of dielectric response 
because it follows quite generally from the /-sum rule jlTj which is a direct consequence 
of particle conservation. At long wavelengths (q — > 0), the plasmon exhausts the /-sum 
rule, and therefore, carries all the spectral weight. In higher dimensions, SPE continua 
start gaining spectral weights as q increases and eventually (q > kp) become the dominant 
elementary excitations. In ID by contrast, the plasmon is the dominant excitation up to 
rather high wave vectors (q > 2kp) and the SPE continua are strongly suppressed. For more 



W(q) 
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(30) 



<9Re[e(g, u))/du 



u>=u) p (q) 
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on S(q, uj) in 3D [17] and 2D we refer to the existing literature [17 , 20 



Fig. 8 shows the calculated weight (Eq. (30)) of the plasmon mode in different dimensions 
as a function of q. In ID (Fig. 8(a)) as q — > W(q) approaches zero because the plasma 
frequency vanishes. It has a maximum strength around q = 2kp for r s = 1.4. For q —>■ oo 
W(q) approaches zero slowly. Note that this absolute vanishing of W(q) for q — > is 
only a manifestation of the vanishing of the ID plasma frequency at long wavelength - the 
plasmon carries all the spectral weight as q — > because the plasmon frequency vanishes in 
the long wavelength limit. In 3D the long wavelength plasmon has a gap, and therefore the 
plasmon weight is finite in 3D even for q = 0. Using Eq. (30), the plasmon weight W(q) 
can be analytically calculated as q — > 0, with the result: W(q) = -Kuo p {q)/2 in the relevant 
dimension. Note that the small weight of the plasmon mode at long wavelengths makes it 
particularly susceptible to collisional damping effects in ID and 2D, as we have discussed 
before. 

In Figs. 9 to 12 we show Re[e(g, uj)], Im[e(g, uj)], and Im[e(g, u;) -1 ] as a function of uj 
for q = 0.2kp. Figs. 9 and 10 show the results for T = (Fig. 9) and T = 0.5Ep (Fig. 
10) with various impurity scattering rates 7, and Figs. 11 and 12 show results for various 
temperatures without any impurity scattering (Fig. 11) and with an impurity scattering rate 
7 = O.lEp (Fig. 12). Both impurity scattering and thermal excitation soften the singularity 
in the dynamical dielectric function at uj + = q 2 + 2q and u;_ = q 2 — 2q (with q measured 
in kp). This causes broadening of the plasmon 5-function peak in the dynamical structure 
factor. As the impurity scattering rate 7 increases the plasmon broadening increases also. 
Finite temperature effects on plasmon broadening are in general small in comparison with 
impurity scattering effects. Even a small impurity scattering (7 m O.lEp) is sufficient to 
cause observable plasmon broadening, but observable thermal broadening of plasmons in 
the clean (7 = 0) system requires a temperature of at least T = 0.5E F . 
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V. HUBBARD APPROXIMATION 



To improve upon RPA, which is what we have so far employed in this paper , and which 



in general is exact only in the long wavelength limit, Hubbard [o,M,22] introduced a simple 
local field correction G(q) to the RPA by approximately summing the ladder exchange dia- 
grams instead of using only the noninteracting polarizability Ho in the dynamical screening: 

/ \ i v c (q)U (q,uj) 

e H {q, uj) = 1- — r, 31 

1 + v c {q)G{q)U (q,uj) 

where 



G{q) = ~ W " , , (32) 
2 v c (q) 

where v c is the Coulomb interaction in the relevant dimension. In the limit q — > 0, G(q) — > 0, 
restoring RPA, which is exact in this limit. Thus in the long wavelength limit the local field 
correction vanishes and Eq. (31) becomes the RPA dielectric function of Eq. (1). The 
simple Hubbard local field corrections in 3D, 2D, and ID can be written |IT| , |2l] , |2"2"|| on the 
basis of Eq. (32) as 

q+kp (33) 



G 3D (g) 



a 



and in ID 



lv e (Jq 2 + kl) 

where v c (q) is now the ID Coulomb interaction in the quantum wire as defined by Eq. (3). 
In Fig. 13 we show the collective plasmon dispersion as calculated by using the Hubbard 
approximation (dashed line) and compared with RPA results (solid line). Dotted lines in Fig. 
13 indicate the boundaries of the electron-hole pair continuum. In the long wavelength limit 
the two approximations match exactly, but as q increases the Hubbard plasmon dispersion 
lies below the RPA plasmon in all dimensions. Thus the Hubbard local field correction 
suppresses the plasmon frequency at finite q in all dimensions. The Hubbard correction, in 
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fact, is larger in lower dimensions. Whether the Hubbard approximation is an improvement 
over RPA in determining the plasmon dispersion is unclear because the Hubbard local field 
correction is purely static whereas the exact local field correction is obviously dynamic, and 
plasmon dispersion is manifestly a dynamical property. We have provided here the plasmon 
Hubbard approximation results only for the sake of completeness. We believe that the 
Hubbard approximation is, in fact, worse than RPA in determining the plasmon dispersion, 
at least in ID and 2D p2| . Within the Tomonaga-Luttinger model Q (i.e. linear bare 
electron energy dispersion and short-range model interaction), RPA is known || to be exact 
because all vertex corrections to the irreducible polarizability vanish ||. For the realistic 
Coulomb interaction and the quadratic electron energy dispersion, which is our ID electron 
gas model in this paper, vertex corrections are not expected to vanish exactly, but should be 
small. Thus, Hubbard approximation, which is a crude estimate of ladder vertex corrections, 
should not be much of an improvement over RPA in ID. Comparison with the experimental 
results of ref. |]IJ] indicates that the experimental ID plasmon dispersion lies above the 
RPA curve whereas we find (Fig. 13) that the Hubbard results lie below the RPA results, 
and therefore inclusion of local field corrections does not by any means improve agreement 
between theory and experiment. In fact, essentially exact quantitative agreement between 
the RPA theory and experiment [|IJ] has recently been obtained || by incorporating a small 
occupancy of a higher ID subband in the calculation. 



VI. CONCLUSION 

We have discussed in this article in considerable details the RPA dynamical response 
of a one dimensional electron gas as confined, for example, in an ultra-narrow semiconduc- 
tor quantum wire structure in the one subband occupancy limit. Special features of our 
work are inclusion of finite impurity scattering and finite temperature effects in the theory 
and a fairly comprehensive discussion of analytic continuation properties of the complex 
dielectric function in the complex frequency plane (necessitated by some confusion on this 
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topic in the recent literature), as well as a detailed comparison between ID plasmon prop- 
erties with corresponding 2D and 3D results. While the motivation for our work has been 
the experimental observation [I] of ID plasmons in GaAs quantum wires via the inelastic 
light scattering spectroscopy, we have not attempted in this paper any comparison between 
theory and experiment because such a comparison already exists in the literature It 
has been established [g| that a complete quantitative and qualitative understanding of all 
aspects of the experimental observations [[[]] requires a two-subband model with the second 
subband occupied slightly by electrons. In the current paper we consider the one-subband 
purely one-dimensional limit which, in fact, explains the basic features of the experimental 
observations |l] very well, as was already noted in ref. Our goal in this paper is to bring 
the understanding of ID dynamical response to the same level as that already existing for 



2D |2],|g and 3D JT^JT^] systems. 

The two most significant features of ID dynamical response are (1) the linear (with small 
logarithmic corrections depending on the details of the confinement potential) wave vector 
dispersion of the long wavelength collective charge density excitation plasmon mode, and (2) 
the total suppression of the electron-hole single particle excitation continua at low energy 
(which arises from the severe ID phase restriction due to energy-momentum conservation). 
The ID plasmon, in fact, exhausts the /-sum rule up to second order in q/kp, and, therefore, 
the ID RPA plasmon dispersion is exact not only in the long wavelength limit (as it is 
trivially in all dimensions) but up to large values of wave vector (q/kp ~ 1). This has earlier 
been pointed out in the literature 0] in a somewhat different context. Our detailed RPA 
treatment of ID dynamical response is thus of considerable quantitative validity because 
RPA turns out to be a more accurate approximation in ID than in 2D and 3D systems. 

Finally, we make some remarks on how the various characteristic peculiarities of ID 
systems, which we neglect uncritically in our work, are expected to affect our RPA theory 
of dynamical response. The three well-known aspects of one dimensionality are Tomonaga- 
Luttinger behavior , Peierls transition and the associated charge density wave instability 
23fl , and Anderson localization |24 |. As mentioned earlier in this article and as has been 



emphasized in ref. 0, a direct implication of the ID Tomonaga-Luttinger behavior is that 
only the bubble diagrams survive || in the irreducible polarizability, making RPA essentially 
an exact approximation for the charge density excitation mode. Thus, what we refer to as 
plasmons throughout this article are the same [[| as the Tomonaga bosons or the charge 
collective excitation modes 0] of a Tomonaga-Luttinger liquid. Much has been written on 
the Peierls transition and associated charge density wave instability in ID systems. We just 
point out that for GaAs quantum wire based on ID electron systems, Peierls transition is 
not a relevant issue as has been emphasized recently in the literature [25|. This is true 



even in the absence of impurity scattering - inclusion of impurity scattering only makes 
Peierls transition and charge density wave instability issue even more irrelevant |H| for one 



dimensional quantum wire electron systems. In ID, any disorder, in principle, completely 



localizes |24| a noninteracting electron system, producing exponentially Anderson localized 
states whereas we have, of course, assumed delocalized plane wave ID electron states in this 
paper. Again, this is not a particularly important issue for our work because in high mobility 



GaAs - based quantum wire systems the calculated localization lengths |2B| are very long 
(> 10 — 100/im), and the system behaves for all practical purpose as a delocalized system. 
Thus the characteristic peculiarities of ID systems, which we have ignored in our work, do 
not in any significant way affect our results or conclusions on the dynamical response of 
quantum wire based ID electron systems. 
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FIGURES 

FIG. 1. (a) The bare bubble diagram, (b) the dressed Green's function, and (c) a typical ladder 
diagram due to impurity scattering. 

FIG. 2. Plasma dispersion in different dimensions. Here the density parameter r s = 1.4 is 
used for all dimensions. The single particle electron-hole continuum excitation regime in ID is 
indicated by the solid curves. 

FIG. 3. The zeros of Re[e(q, u — ia)\ and lm[e(q,u — ia)\ in the complex uj — a plane, where 
a indicates the damping of the plasmon mode, for q = 0.3/cf- Point P corresponds to the location 
of the plasmon: Re[e] = Im[e] = 0. (a) Clean system with the impurity scattering rate 7 = 0. (b) 
Dirty system with impurity scattering rate 7 = O.lEp- Note that P lies on the a = line in (a), 
but is shifted to finite a in (b). 

FIG. 4. ID plasmon dispersion is shown for various impurity scattering rates 7 as given in the 
figure. Here uj > gives the dispersion of the mode and u < the damping of the mode. The 
overdamping of the plasmon for q < q c is clearly seen. Parameters are a = lOOA and r s = 1.4 
which corresponds to the ID density n = 0.56 x 10 6 cm _1 for GaAs. 

FIG. 5. Plasmon dispersion and damping in (a) 2D and (b) 3D (r s = 1.4 in both). 

FIG. 6. (a) ID plasmon dispersion and damping without any impurity scattering for various 
temperatures as shown, (b) ID plasmon dispersion and damping at a finite temperature T = O.bEp 
and for various impurity scattering rates 7 as shown (r s = 1.4). 

FIG. 7. Im[e(g, lo)^ 1 ] in ID as a function of frequency u for different wave vectors q. The 
(5-function peaks are the plasmon modes with their spectral weights written above the peaks. 
The strength of the electron-hole continuum (the broad incoherent bump in each figure) has been 
enhanced by a factor ten for visual clarity. 
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FIG. 8. Calculated weights W(q) of the plasmon mode in different dimensions: (a) ID, (b) 2D, 
and (c) 3D. 

FIG. 9. Shows for a fixed wave vector q = 0.2/ci?: (a) Re[e(q,uj)]; (b) Im[e(g, oj)]; and (c) 
Im[e(g, a;) -1 ] as a function of frequency at T = with various impurity scattering rates 7. In (c) 
the electron-hole pair continuum for 7 = (solid line) has been enhanced by a factor 10 for the 
sake of clarity. The vertical arrow in (c) denotes the plasmon 5-function peak. 

FIG. 10. Shows for a fixed wave vector q = 0.2/c^ the calculated (a) Re[e(q, 00)]; (b) Im[e(g,cj)]; 
and (c) Im[e((7, oj)^ 1 ] as a function of frequency to at a finite temperature T = 0.5Ep with various 
impurity scattering rates 7 as shown. 

FIG. 11. Shows for a fixed wave vector q = 0.2/c^ the calculated (a) Re[e(g, oj)]; (b) lm[e(q, oj)]; 
and (c) lm[e(q,uj)~ 1 } as a function of frequency uj without any impurity scattering (7 = 0) for 
various temperatures as shown. 

FIG. 12. Shows for a wave vector q = 0.2/c^p the calculated (a) Re[e(q, uj)]; (b) Im[e(g,cj)]; 
and (c) Im[e(</, oj)" 1 ] as a function of frequency to with an impurity scattering rate 7 = O.lEp for 
various temperatures as shown. 

FIG. 13. Collective plasmon dispersion within the Hubbard approximation (dashed lines) and 
within RPA (solid lines) for different dimensions: (a) ID; (b) 2D; and (c) 3D. Dotted lines indicate 
the boundaries of the electron-hole pair continua. 
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